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' Abstract 

F ^ \ The Painleve test is very useful to construct not only the Laurent series solutions of 

|"T i ■ systems of nonlinear ordinary differential equations but also the elliptic and trigono- 

| metric ones. The standard methods for constructing the elliptic solutions consist of 

■ two independent steps: transformation of a nonlinear polynomial differential equation 

into a nonlinear algebraic system and a search for solutions of the obtained system. It 
has been demonstrated by the example of the generalized Henon-Heiles system that 
the use of the Laurent series solutions of the initial differential equation assists to solve 
the obtained algebraic system. This procedure has been automatized and generalized 
on some type of multivalued solutions. To find solutions of the initial differential equa- 
' tion in the form of the Laurent or Puiseux series we use the Painleve test. This test 

can also assist to solve the inverse problem: to find the form of a polynomial potential, 
which corresponds to the required type of solutions. We consider the five-dimensional 
gravitational model with a scalar field to demonstrate this. 
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1 Introduction 



The investigations of exact special solutions of nonintegrable systems play an important 
role in the study of nonlinear physical phenomena. There are a few methods to construct 
solutions in terms of rational, hyperbolic, trigonometric or elliptic functions (all such 
^ . functions are solutions of the first order polynomial differential equations) [21 El El E3 El 

E3E31E3I33IBHI33I3B1I13I1H- These methods (at least, some of them) use information 
about the dominant behavior of the initial system solutions in the neighbourhood of their 
singular points, but do not use the Laurent series representations of them. In [201 a new 
method has been developed as an alternative way to construct elliptic and elementary 
solutions and it has been shown that the use of the Laurent series solutions allows to 
extend the set of solutions ^U] or, on the contrary, to prove non-existence of elliptic 
solutions |21| I42j . The Laurent series solutions give the information about the global 
behavior of the differential system and assist in looking not only for global solutions but 
also for first integrals |24j . 

The direct algebraic method is the use of the first order polynomial differential equation, 
which solutions are elementary or elliptic functions, to transform the initial differential 
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system in a nonlinear algebraic system in coefficients of this first order equation and 
parameters of the initial system. It has been proved that any finite system of algebraic 
equations can be solved by the Grobner (Groebner) basis method ^2] 5 but calculations for 
a sufficiently complex algebraic system can be very difficult and expensive. In this paper 
we show that the use of the Laurent series solutions gives additional algebraic equations 
and allows calculations to be simplified. These equations are linear in coefficients of 
the first order equation and nonlinear, maybe even nonpolynomial, in parameters of the 
initial system. Note that one maybe should fix some of these parameters to construct 
the Laurent series solutions. Therefore, in contrast to the Grobner basis method, the 
additional equations may be not consequences of the initial algebraic system. 

To manifest the main properties of this method we consider the generalized Henon- 
Heiles system, for which analytic and Laurent series solutions have been found in |35| I36| 

EH]. 

2 Painleve analysis 

We interpret solutions of a system of ordinary differential equations (ODE's) as analytic 
functions, maybe with isolated singular points. A singular point of a solution is said 
critical (as opposed to noncritical) if the solution is multivalued (single-valued) in its 
neighbourhood and movable if its location depends on initial conditions A singular 
point is said algebraic if the solution has a finite number of values in its neighbourhood. 
In the opposite case a singular point is either transcendental or essentially singular one. 

The general solution of an ODE of order N is the set of all solutions mentioned in 
the existence theorem of Cauchy, i.e. determined by the initial values. It depends on TV 
arbitrary independent constants. A special solution is any solution obtained from the 
general solution by giving values to the arbitrary constants. A system of ODE's has the 
Painleve property if its general solution has no movable critical singular point |31j . 

The Painleve test is any algorithm, which checks some necessary conditions for a dif- 
ferential equation to have the Painleve property. The original algorithm, developed by 
P. Painleve and used by him to find all the second order ODE's with the Painleve prop- 
erty, is known as the a- method. The method of S.V. Kovalevskaya |2f)] is not as general 
as the a-method, but much more simple. In 1980, developing the Kovalevskaya method 
further, M.J. Ablowitz, A. Ramani and H. Segur ,1 constructed a new algorithm of the 
Painleve test for ODE's. This algorithm appears very useful to find solutions as a formal 
Laurent series |39j . First of all, it allows to determine the dominant behavior of a solution 
in the neighbourhood of the singular point to- If the solution tends to infinity as (t — to)! 3 ', 
where (3 is a negative integer number, then substituting the Laurent series expansion one 
can transform a nonlinear differential equation into a system of linear algebraic equations 
on coefficients of the Laurent series. All solutions of a system depend on the parameter 
to, which characterizes the singular point location. If a single-valued solution depends on 
other parameters, then some coefficients of its Laurent series have to be arbitrary and 
the matrices in the corresponding linear systems should have zero determinants. The 
characteristic numbers associated with such systems (named resonances or Kovalevskaya 
exponents) are thus determined using the Painleve test. In such a way we obtain solutions 
only as formal series, but really we will use only a finite number of these series coefficients, 
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so, we do not need the convergence of these series. At the same time for some noninte- 
grable systems, in particular, the generalized Henon-Heiles system, the convergence of the 
Laurent- and psi-series solutions has been proved |28| . 

Classical results allow one to construct the suitable form of the first order autonomous 
polynomial equation. A theorem due to Painleve |31[ I14j proves that solutions of the 
equation 

m 

J2Pk(y(t),t)y t (t) k = o, (1) 

k=Q 

where y t (t) = dy Iy , m is a positive integer number and Pk are polynomials in y(t) and 
analytic functions in t, have neither movable transcendental nor movable essential singular 
point. A theorem of Fuchs |181 114j shows that if eq. Q has no critical movable singular 
point, then the degree of Pk as a polynomial in y(t) is no more than 2m — 2k, in particular, 
P m does not depend on y. 

Therefore, the necessary form of a polynomial autonomous first order ODE with the 
single- valued general solution is 

m 2m— 2k 

E E h jkV ] yt=o, h 0m = i, (2) 

k=0 j=0 

in which m is a positive integer number and hjk are constants. A theorem of Briot and 
Bouquet [SJ proves that the general solution of © is either an elliptic function, or a 
rational function of e iX , 7 being some constant, or a rational function of x. Note that the 
third case is a degeneracy of the second one, which in turn is a degenerate case of the first 
one. 



3 The Henon— Heiles system 

To compare methods for the construction of elliptic solutions let us consider the generalized 
Henon-Heiles system with an additional non-polynomial term, which is described by the 
Hamiltonian: 

H = ~(x t 2 + y 2 t + Aire 2 + A 2 y 2 ) + x 2 y - - y 3 + ^ (3) 
and the corresponding system of equations of motion: 

x t t = - \ix -2xy + A,, 

x 6 (4) 

Vu = - \2y-x 2 + Cy 2 , 

where subscripts denote derivatives: x tt = and y tt = 4r|f ; Ax, A 2 , \i and C are arbitrary 
numerical parameters. Note that if A 2 7^ 0, then one can put A 2 = si<7n(A 2 ) without the 
loss of generality. 

If C = 1, Ai = 1, A 2 = 1 and [i = 0, then Q is the original Henon-Heiles system |19j . 
Emphasizing that their choice of potential does not proceed from experimental data, Henon 
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and Heiles have proposed this system to study the motion of a star in an axial-symmetric 
and time-independent potential, because, on the one hand, it is analytically simple and 
this makes the numerical computations of trajectories easy, but, on the other hand, it 
is sufficiently complicated to give trajectories which are far from trivial. Subsequent 
numerical investigations E| show that in the complex t-plane singular points of solutions 
of the Henon-Heiles system group in self-similar spirals. It turns out that there are 
extremely complicated distributions of singularities, forming a boundary, across which the 
solutions can not be analytically continued. 

By the use of Painleve analysis the following integrable cases of the generalized Henon- 
Heiles system have been found: 

(i) C = -1, Ai = A 2 , 

(ii) C = —6, Ai, A2 arbitrary, 

(iii) C = -16, Ai = A2/I6. 

In all above-mentioned cases system (J1J) is integrable for any value of [i. The general 
solutions of the Henon-Heiles system are known only in integrable cases |37[ I38( II lj , in 
other cases not only four-, but even three-parameter exact solutions have yet to be found. 
The generalized Henon-Heiles system has attracted enormous attention over the years and 
has used as a model in astronomy j3U] and in gravitation [251 EH] • 



4 Methods for construction of elliptic solutions 

4.1 The construction of the nonlinear algebraic system 

The function y, solution of system (fijl. satisfies the following fourth-order polynomial 
equation: 

y mt = (2C-8)y tt y-(4X 1 + X 2 )y tt + 2(C + l)y 2 + 

2QC ( 5 ) 
+ ^-y 3 + (4CAi - 6A 2 )y 2 - 4A!A 2 y - AH. 

We note that the energy of the system H is not an arbitrary parameter, but a function of 
initial data: y , y ot , y ott and y ota - The form of this function depends on /i: 

rr vlt + Vl C vl , A i + 2 Vo » \ 1 (A 2 y 0t + 2Cy y 0t - y mt ) 2 + \x 

This formula is correct only if x = Cy\ — A 2 y — y m 7^ 0. If x = 0, which is possible 
only for /u = 0, then we can not express x ot through y , y ot , y ott and yQ ttt , so H is not 
a function of the initial data. If y ottt = 1Cy§y§ t — A 2 y 04 , then eq. © with an arbitrary 
H corresponds to system (jlj) with /j, = 0, in opposite case eq. (JSJ) does not correspond to 
system 

To find a special solution of eq. ® one can assume that y satisfies some more simple 
equation. For example, there exist solutions in terms of the Weierstrass elliptic function (or 
a degenerate elliptic function), which satisfy the following first-order differential equation: 



y\ = Ay 3 + By 2 + Cy + D, 



(6) 
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where A, B, C and D are constants to be determined. 
E.I. Timoshkova [35] generalized eq. ©: 

y 2 t = A 4 y 3 + A 3 y 5 ' 2 + A 2 y 2 + A lV 3 / 2 + A y + D, 



(7) 



(Aj are constants) and found new one-parameter solutions of the Henon-Heiles system in 
two nonintegrable cases (C = —4/3 or C = — 16/5). The Timoshkova's method allows 
one to seek both single- and multivalued solutions. For the generalized Henon-Heiles 
system new solutions (i.e. solutions with ^3 7^ or A\ 7^ 0) are derived only at D = 0. 
These solutions are single- valued, because the substitution y(t) = g(t) 2 gives: 



Qt = t (Mq a + A 3 g 3 + A 2 q 2 + A ie + A)) • 
Following we use the substitution 
y(t) = g{t) 2 + P , 
where Pq is a constant, and transform eq. © into 

&ttt0 = -4 % ^ t -3^ + 2(C-4) % ^ + (2P (C-4)-4A 1 -A 2 )^^ + 



(8) 



0) 



+ 2(3C - 2) Q 2 te 2 + (2CP - 4Ai - 8P - A 2 ) ft 2 + ^-Cg 6 + 
+ (2CA a + WCP - 3X 2 )g 4 + 2(2X 1 CP + 5CP 2 - A X A 2 - 3P A 2 )£ 2 + 

+ -z;CF$ + 2AiCP 2 - 3P 2 A 2 - 2A!A 2 P - 2H. 

If g(t) satisfies (jSJ), then eq. (jTUI) is equivalent to the following system 

{3A A + 4) (2C - 3A4) = 0, 
A 3 (9C-21At - 16) = 0, 

96A A CP - 240A A A 2 - 192AjAi - 384vl4Po - 48A*A 2 - 

- 105Aj + 128A 2 C - 192A 2 + 128CAi + 640CP - 192A 2 = 0, 

40A 3 CP - 90^4^1 - 65^ 3 ^ 2 - 8OA3A1 - 

- 160A 3 P - 2(L4 3 A 2 + 56CAi - 64Ai = 0, 

16^ 2 CP - S6A4A0 - 21.43.4i - 8A 2 2 - 32A 2 X 1 - 64A 2 P - 8X 2 A 2 + 
+ 24CA) + 64AiCP + 160CP 2 - 16A - 32A X A 2 - 96P A 2 = 0, 

IOA3A0 + (5A 2 + 8CP - I6A1 - 32P - 4A 2 )^i = 

and the equation for the energy H: 
1 



10 



(10) 



(11) 



H 



384 



96CA P - 48A A 2 + 384CAiP,f + 640CP d - 9 A 



192A Ai - 384A P Q - 48A A 2 - 384AiA 2 P - 576A 2 P ' 



(12) 
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It has been proposed J3| to seek solutions as polynomials with three arbitrary coeffi- 
cients 

y(t) = P 2 g(t) 2 + P 10 (t) + P , (13) 
where the function g(t) satisfies eq. (JHJ). The function 



^(*»-f)- (14> 



satisfies eq. (jSJ) as well, therefore, one and the same function y(t) corresponds to a two- 
parameter set of coefficients A^ [i = 0..4) and P& (k = 0..2), so we always can put 

P 2 = 1 and Pi = 0, (15) 

without loss of generality. 

In the general case, when one uses Fan's technique ^7] (see, for example, [33]) with 
polynomials of degree L in solutions of eq. (jHJ) or eq. @, it is possible to simplify calcu- 
lations, putting 

P L = 1 and P L -i = 0. (16) 

In particular, the consideration of linear combinations P\y + Po instead of the function y 
is useless. 

System tf 1 If) has been solved by REDUCE using the standard function solve and the 
Grobner basis method jHS]. The goal of this paper is to show that the use of the Laurent 
series solutions allows to obtain some solutions of (|11|) solving only linear systems and 
nonlinear equations in one variable. 

We cannot use this method for arbitrary C, because the Laurent series solutions are 
different for different C. So, first of all we have to fix value of C. Note that if ^3 = A\ = 
then y satisfies eq. (jHJ) with D = 0, hence, such solutions are known. The remaining 
problem can be separated on two cases: A3 7^ or ^3 = and A\ 7^ 0. If ^3 7^ 0, then 
from two first equations of system 1)11(1 it follows that 

C = -* and * = or O = - ™ and A 4 — — — . 

3 3 5 15 

Surely, if some solution with ^3 = corresponds to the above-mentioned values of C, then 

it can be found as well. 

4.2 Construction of linear algebraic system 

Let choose C = — 4/3 (the case C = — 16/5 one can consider analogously). After 
substitution C = — 4/3 and A4 = — 4/3 in (|TTj) we obtain 

512P + 128vl2 - 256Ai + 384A 2 + 315^ = 0, 

56^i + (640P + 195A 2 + 240Ai + 60X 2 )A 3 = 0, 

63^3^1 + 24 (a 2 + 4A1 + A 2 + y P ) A 2 + (17) 

+ 256AiP + 640P 2 + 96AiA 2 + 288P A 2 = 0, 

/ 128 \ 

10A 3 A Q +UA 2 - — P - I6A1 - 4A 2 J A x = 0. 
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If we consider system (|T7|) separately from the differential equations © and (|TU|) , from 
which it has been obtained, then it would be difficult to solve system (|17p without the 
use of the Grobner basis method. At the same time from equations (jHJ) and we can 
obtain additional information, which assists us to solve system (|17|) . 

Let us construct the Laurent series solutions for eq. I)10|l . The method of construction 
of the Laurent series solutions for the generalized Henon-Heiles system has been described 
in detail in [22] • For eq. (jlOj) with C = — 4/3 we obtain that solutions have singularities 
proportional to 1/t and the values of resonances are — 1 (corresponding to the arbitrary 
parameter to), 1, 4 and 10. The Laurent series solutions are (we put to = 0): 

^ ,V ^ + co + ^(3A 2 -2Ai+4P + 62cg)t + ... ) , (18) 



t 24 



where 



±a/161700Ai - 121275A 2 ± VH55(5481A| - 12768AiA 2 + 8512A^) 

C0 = 2311 • (19) 

Two signs " ±" in (|19|) are independent. At the same time, functions g and — g correspond 
to one and the same function y, so there are four different Laurent series solutions. The 
coefficients C3 and eg are arbitrary. To find any number of coefficients, we should solve 
only two nonlinear equations in one variable and linear equations in one variable. 

The algorithm of the construction of elliptic solutions from the Laurent series solutions 
is the following 29 : 

• Choose a positive integer m and define the first order ODE (j2j) with unknown coef- 
ficients hjk- 

• Compute coefficients of the Laurent series g. The number of coefficients has to be 
greater than the number of unknowns. 

• Substituting the obtained coefficients, transform eq. © into a linear and over- 
determined system in hjk with coefficients depending on arbitrary parameters. 

• Exclude hjk and obtain the nonlinear system in parameters. 

• Solve the obtained system. 

To obtain the explicit form of the elliptic function, which satisfies the known first order 
ODE, one can use the classical method due to Poincare, which has been implemented in 
Maple as the package "algcurves" |2*U] . 

On the first step we choose eq. (j2J), which coincides with eq. (jHJ). It means that m = 2, 
all hj\ are equal to zero and all hjo = — Aj/4. After the second and the third steps 
we obtain a linear system in Aj. The package of computer algebra procedures, which 
transforms the first order equation into a such system of algebraic equations has been 
written |4f)j in Maple. 

The obtained system has the triangular form and is linear in H, C3 and eg as well. From 
the first equation we obtain anew that A4 = — 4/3. From the second equation it follows: 

^3 = ^c , (20) 
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and so on: 

„2 



A 2 = - 70c z - 3A 2 + 2Ai - 4P , (21) 



.1, = ( yP - 60Ai + 50A 2 + 1300c 2 , ) r n . (22) 



(23) 



, r 21535 4 /565, 405 , 245 „ \ 2 
4, = - 40z\/3c3 - ^-cg + ( — Ai-— A2--3- P0UI + 

7 21 2 7 2 7 7 7 2 

+ — Ai A 2 — Ao — Ai + — Ai Pn — — AoPn — —Pn ■ 

4 16 2 12 1 3 2 u 3 

From the next equation of the system we obtain C3 and, finally, 

15645 4 /1545 x 1495 \ 2 537 x2 

4 = — A 2 -465P -— ^Jc^— Af - 

- ^ Al A 2 + ^A| + 19^0-^0-^. 
20 80 2 2 u 3 

Substituting the values of Af., which correspond to one of the possible values of cn, in 
system (|17[) we obtain that it is satisfied for all values of Ai, A 2 and Pq, so we do not need 
to solve the nonlinear equations. Therefore we settle the nonlinear algebraic system ()17|) . 
solving only linear equations and nonlinear equation in one variable. We have used the 
values of only six coefficients of the Laurent series solutions. Note that, for cq and — c$ we 
obtain one and the same values of A±, Ai and An., and the opposite values of ^3 and A\. 
From (JHJ) it follows that these solutions correspond to ±g(t), and, hence, give one and the 
same function y(t). Therefore not four, but two different elliptic (or degenerated elliptic) 
solutions of eq. © have been found. The elliptic solutions y(t) are fourth-order elliptic 
functions and can be expressed in terms of the Weierstrass elliptic function p(t — to): 

where constants a, b, c, d and periods of p(t) are determined by A^. The parameters to 
and Pn, which defines the energy of the system, are arbitrary. Solutions of this type exist 
in both above-mentioned nonintegrable cases: C = — 16/5 and C = — 4/3. The full list 
of solutions is given in |36j , 

5 Multivalued solutions 

From the Painleve theorem it follows that a solution of an autonomous polynomial first or- 
der differential equation can have only two types of singular points: poles and algebraic 
branch points. So, the most general type of formal series solutions is the Puiseux series. 
Let us generalize the above method by searching for solutions, which can be expanded in 
Puiseux series, as follows: 



'k/ q t k/q - (26) 
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We seek solutions polynomial 

y = P L p L + P L - lP L ^ + P L -2P L ~ 2 + Pl-3P L ~ 3 + • • • + Po- (27) 
From Q26JI it follows that 

CO 

P = E C i/^ /q - ( 28 ) 
i=-i 

We assume that p(t) satisfies the following equation: 

m (g+l)(m— k) 

E E h jk p>p t k = 0, h 0m = l, (29) 

k=0 j=0 

where q is a natural number. 

To simplify calculations one can put Pl = 1 and Pl~\ = without loss of generality, 
because 

^(p-^j/lfa (30) 

satisfies ()29|) as well. The construction of the corresponding algebraic system has been 
automatized in Maple jlj. We plan to automatize similar construction for an arbitrary 
autonomous polynomial first order equation. 

6 Five— dimensional gravitational model with a scalar field 

To show how the analysis of singular behavior of solutions can assist to find the form of 
a potential, let us consider a five-dimensional model of a gravitational field ^3E]> w ^h 
the action given by 



S = J d 4 xdr^J\ det g^ u \ I j - d^(pd u (p - V((p) 



(31) 



where M is the full five-dimensional space-time. The most general metric with four- 
dimensional Poincare symmetry is 

ds 2 = e 2A{r) ( - dx\ + dxj + dx\ + dxfj + dr 2 . (32) 

Assuming the scalar field depends only on r, so that <p = 4>{ r )i the independent equations 
of motion are 

' J' = 

1-1 ( 33 ) 

J 2 = \ W) - \v{4>), 

where J = A' = $ and 0' = 
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Similar systems of equations of motion arise also in cosmological models with a spatially 
flat Friedmann metric and a scalar field :4 or a phantom scalar field [3]. 

Let us analyse the correspondence between 0(r) and the potential V(4>). If at a singular 
point (fixed at the origin, without loss of generality) 

Mr) ~ — , 
then from (|33|) it follows that 

J'~^To J~^T—T =>• ~ V f — ] ~ -7^X9. (34) 

It means that solutions with simple poles can be obtained only if the degree of the 
polynomial potential V(4>) is equal to six. For example, let <p{r) = tanh(r), for real r this 
function has no singular point, but system (|3*3*)) is autonomous, so if 4>{r) is a solution, 
then 4>{r — ro), where ro is an arbitrary complex constant, is a solution as well. Then there 
exists ro such that tanh(r — ro) tends to infinity as r — ► 0, hence, the function tanh(r) 
cannot be a solution for any fourth degree polynomial potential. In ^3] the explicit form 
of the sixth degree polynomial potential V(cp), which corresponds to tanh(r) has been 
found. Analogously one can show that if solutions tend to infinity as 1/r 2 , then the degree 
of V{cj)) is equal to 5. If solutions tend to infinity as l/r fc , where A; is a natural number 
greater than two, then from Q34|) we obtain that V(4>) cannot be a polynomial. 

The polynomial potentials, which degrees are more than 6, correspond to multivalued 
solutions. If solutions tend to infinity as r _2//fc , then from (|34|) it is follows that the degree 
of V(4>) has to be equal to 4 + k. For example, if a solution in the neighbourhood of its 
singular points tends to infinity as 1/y/r, then k = 4 and the degree of the corresponding 
potential V(4>) is equal to 8. 

In conclusion of this section we say a few words about explicit solutions. Following ^3] 
we assume that J(r) is a function of 4>: 

J(r) = ~ \w(<t>(r)). (35) 



It is straightforward to verify that system (|33|) is equivalent to 

d(j)(r) _ 1 dW(<j>) 
dr ~ 2 d6 ' 



(36) 



^) 2 -~^) 2 -^) = 0. (37) 

Unfortunately the nonautonomous equation (|37|) can not be solved analytically. For 
polynomial V(4>) it may be possible to find only special solutions, e.g. W(4>) in polynomial 
form. For example, W{4>) cannot be a polynomial if V(<p) = ((j) 2 — l) and we don't know 
a solution for this potential. At the same time, if the potential V{4>) = A<j) 4 + Bcp 2 + C 
has no double roots, then eq. (|37|) has a polynomial solution. 
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Contrary to a scalar field theory without gravitational field, there is not a one-to-one 
correspondence between the form of the scalar field cj)(r — tq) and potential V((f>). The 
form of the scalar field is defined by so one can add a constant to W(<p) and obtain 
new V{4>) for the same 4>(r). On the other hand, for given V((p) we have not one-, but 
two-parameter set of functions 4>{r). 



7 Conclusion 

The Laurent series solutions are useful to find elliptic or elementary solutions in the ana- 
lytic form. The method, proposed in has been automatized and generalized on some 
type of multivalued solutions. It converts the local information into the global one and 
can be used not only as an alternative to the standard method, but also as an addition to 
it, which assists to find solutions of the obtained algebraic system. We have demonstrated 
that one can find elliptic solutions of the generalized Henon-Heiles system solving only 
linear equations and nonlinear equations in one variable, instead of nonlinear system 
At the same time, to use this method one has to know not only an algebraic system, but 
also the differential equations from which this system has been obtained. To find solutions 
of the initial ODE in the form of the Laurent or Puiseux series we use the Painleve test. 
This test can also assist to solve the inverse problem: to find the degree of a polynomial 
potential, which corresponds to the required type of solutions. 
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